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Abstract 

The effect of event background fluctuations on charged particle jet reconstruction in Pb-Pb colli- 
sions at = 2.76 TeV has been measured with the ALICE experiment. The main sources of non- 
statisticai fluctuations are characterized based purely on experimental data with an unbiased method, 
as well as by using single high p t particles and simulated jets embedded into real Pb-Pb events and 
reconstructed with the anti-fe t jet finder. The influence of a low transverse momentum cut-off on parti- 
cles used in the jet reconstruction is quantified by varying the minimum track p t between 0.15 GeV/c 
and 2GeV/c. For embedded jets reconstructed from charged particles with p t > 0.15GeV/c, the 
uncertainty in the reconstructed jet transverse momentum due to the heavy-ion background is mea- 
sured to be 11 .3 GeV/c (standard deviation) for the 10% most central Pb-Pb collisions, slightly larger 
than the value of ll.OGeV/c measured using the unbiased method. For a higher particle transverse 
momentum threshold of 2 GeV/c, which will generate a stronger bias towards hard fragmentation in 
the jet finding process, the standard deviation of the fluctuations in the reconstructed jet transverse 
momentum is reduced to 4.8-5.0 GeV/c for the 10% most central events. A non-Gaussian tail of 
the momentum uncertainty is observed and its impact on the reconstructed jet spectrum is evaluated 
for varying particle momentum thresholds, by folding the measured fluctuations with steeply falling 
spectra. 



*See Appendix A for the list of collaboration members 
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1 Introduction 

High energy heavy-ion collisions explore strongly interacting matter under extreme conditions of energy 
density, where lattice QCD predicts a phase transition to a new state of matter above a critical value of 
about 1 GeV/fm (H. In this new state, called the Quark-Gluon Plasma (QGP), quarks and gluons rather 
than hadrons are expected to be the dominant degrees of freedom over length scales larger than that of a 
nucleon. Experiments studying the collision of heavy nuclei at high energy at both the Relativistic Heavy 
Ion Collider (RHIC) E [2 13 S, and recently at the Large Hadron Collider (LHC) IE 13 gill, have made 
several key observations that point to the formation of a hot, dense and strongly coupled system, possibly 
the QGP 

Hard (large momentum transfer Q 2 ) probes are well calibrated tools to study the properties of the matter 
created in such collisions. The scattered partons generated in a hard momentum exchange are created in 
the initial stages of the heavy-ion collision, with production rates that are calculable using perturbative 
QCD, which can be compared to the same measurements in proton-proton collisions. The scattered 
partons then propagate through the medium, where their fragmentation into observed jets of hadrons is 
expected to be modified relative to the vacuum case by interactions with the medium (jet quenching) 
ifTOlfTTI . This modification of parton fragmentation provides sensitive observables to study properties of 
the created matter. 

Jet quenching has been observed at RHIC |[T2l [131 [I4l [TBI and at the LHC lfT6l via the measurement of 
high p t hadron inclusive production and correlations, which are observed to be strongly suppressed in 
central A-A collisions compared to a scaled pp reference. These high p x hadron observables have been 
the major tool for measuring the energy loss of hard scattered partons and thereby the properties of the 
medium, but they provide only indirect and biased information on the parton evolution in the medium. 
The aim of full jet reconstruction is to measure jet modifications due to energy loss in an unbiased way 
lfT71[T8l . Already first measurements of reconstructed jets in heavy-ion collisions at the LHC showed an 
energy imbalance between back-to-back dijets, which is attributed to jet quenching |9l[T9ll. 

Jet reconstruction in the complex environment of a heavy-ion collision requires a quantitative under- 
standing of background-induced fluctuations of the measured jet signal and the effects of the underlying 
heavy-ion event on the jet finding process itself. Here, region-to-region background fluctuations are the 
main source of jet energy or momentum uncertainty and can have a large impact on jet structure ob- 
servables, such as the fraction of energy inside the jet core or the shape of the jet, and will distort the 
measured jet energy balance even in the absence of medium effects lf20l . 

In this paper the measurement of jet transverse momentum fluctuations due to the background in heavy- 
ion collisions is reported and its sources are identified, based on jet reconstruction using charged particles 
with varying minimum track p t . For this purpose three methods are employed to probe the measured 
Pb-Pb events: fixed area (rigid) cones placed randomly in the acceptance, the simulation of high-p t 
single tracks or full jets from pp collisions. Rigid cones enable the identification of contributions to the 
fluctuations in an unbiased fashion, while single tracks and embedded jets explore the interplay between 
the jet finding process, the underlying event, and the jet fragmentation pattern. 

2 Detector Description and Track Selection 

The data presented here were collected by the ALICE experiment ll2~TTl in the first Pb-Pb run of the 
LHC in November 2010, at a collision energy of y/s^ = 2.76 TeV. This analysis is based on minimum- 
bias events, triggered by two forward VZERO counters and the Silicon Pixel Detector (SPD) ll22l . A 
description of the minimum-bias trigger can be found in |6). The VZERO trigger counters are forward 
scintillator detectors covering a pseudo-rapidity range of 2.8 < r\ < 5.1 (VOA) and —3.7 < r\ < —1.7 
(VOC). The sum of VZERO amplitudes is also used as a measure of event centrality ll23l . The SPD 
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consists of two silicon pixel layers at a radial distance to the beam line of r = 3.9 cm and r = 7.6 cm. 

To ensure a uniform track acceptance in pseudo-rapidity r\ , only events whose primary vertex lies within 
±8 cm from the center of the detector along the beam line are used, resulting in 13.3 M minimum-bias 
Pb-Pb events for this analysis. 

Charged particle tracking is carried out using the Time Projection Chamber (TPC) B4l and the In- 
ner Tracking System (ITS) [22], located in the central barrel of the ALICE experiment within a 0.5 T 
solenoidal magnetic field and covering the full azimuth within pseudo-rapidity |tj| < 0.9. The ITS con- 
sists of six cylindrical layers of silicon detectors, with distances from the beam-axis between r = 3.9 cm 
and r = 43 cm. The ITS layers measure track points close to the primary vertex, with the two innermost 
layers (SPD) providing a precise measurement of the primary vertex position. The TPC, a cylindrical 
drift detector surrounding the ITS, is the main tracking detector in ALICE. The TPC inner radius is 85 cm 
and the outer radius is 247 cm, with longitudinal coverage —250 < z < 250 cm. It provides a uniformly 
high tracking efficiency for charged particles. The high precision of the ITS and the large radial lever 
arm of the TPC provide a good momentum resolution for combined (global) tracks. 

For this dataset the ITS has significantly non-uniform efficiency as a function of azimuthal angle </> and 
pseudo-rapidity v\ . In order to obtain high and uniform tracking efficiency together with good momentum 
resolution, two different track populations are utilized: (i) tracks containing at least one space -point 
reconstructed in one of the two innermost layers of the ITS (78% of all accepted tracks), and (ii) accepted 
tracks that lack the position information close to the beam-line. Here, the primary vertex is added to the fit 
of the track which modifies the reconstructed curvature of the charged track in the magnetic field. Since 
the majority of the tracks originates from the primary vertex this constraint improves the momentum 
resolution. Both track types have transverse momentum resolution of a(p t )/p t 1% at 1 GeV/c. For 
the majority of tracks the resolution at p t = 50GeV/c is o(p t )/pt ~ 10%, only tracks having fewer 
than three reconstructed space points in the ITS (about 6% of the total population) have a resolution at 
50GeV/c of o(pt)/pt ~ 20%. 

Tracks are accepted for p t > 0.15GeV/c and |t]| < 0.9. The tracking efficiency at p t = 0.15GeV/c is 
50%, increasing to 90% at 1 GeV/c and above. Tracks with measured p t > lOOGeV/c are accepted at 
the tracking stage, but jets containing them are rejected from the analysis to reduce the influence of fake 
tracks and limited tracking resolution at very high p t . 

3 Jet Reconstruction and Background Subtraction 

Charged particle jet reconstruction and estimation of the background employ the sequential recombina- 
tion algorithms anti-& t and k t from the FastJet package ll25l . The clustering starts with the list of tracks 
that satisfy the quality, acceptance, and /? t -cuts, with no pre-clustering or grouping of tracks. A list of jet 
candidates (anti-& t ) or clusters (k t ) is generated, with direction and transverse momentum given by the 
/7 t -weighted average of (t],0) of the individual constituents and the scalar sum of their p u respectively. 
The distance parameter that determines the terminating condition for the clustering is chosen as R = 0.4, 
which is a common value for reconstruction of jets in heavy-ion collisions IT8l l9l l2"6ll . As proposed in 
E71 . the clusters found by the k t algorithm are used to estimate the event-wise background /? t -density per 
unit area, p, defined as the median value of the ratio p{ ec /A iec for all considered fc t -clusters. A rec is the 
area of the reconstructed cluster in the (17, 0)-plane calculated by the active ghost area method of FastJet 
ll28l . with a ghost area of 0.005. To minimize the influence of the track acceptance interval on p, only 
reconstructed clusters with |t]| < 0.5 have been used. In addition, the two clusters with the largest p[ ec 
(leading) in the full acceptance of |t]| < 0.9 are excluded from the calculation of the median to further 
reduce the influence of true jets on the background estimate [29 ]. The jet population reconstructed by the 
anti-^t algorithm is used as signal jets. Their p t is corrected for the background /? t -density in each event 
using the jet area A-j et ' rec with pf = p{ eUec — p .AJ et ' rec . Signal jets are only considered for \t] \ < 0.5. 
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Figure 1 : Dependence of charged particle background p t density p on uncorrected multiplicity of tracks 
used for jet finding (|t]| < 0.9). The dotted line is a linear fit to the centroids in each multiplicity bin. 
The insets show the projected distributions of p and raw multiplicity for the 10% most central events. 
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Table 1 : Average and standard deviation of the event- wise charged particle p t density p for three choices 
of minimum particle p t and two centrality bins. The quoted uncertainties are purely statistical. 



The average transverse momentum of tracks (p t ) and the total charged multiplicity are global observables 
that are closely related to the value of p, though the determination of p uses varying phase-space intervals 
(with typical areas in the (T] , 0)-plane of 7iR 2 ) and suppresses hard jet contributions by using the median 
of the distribution. Figure [T] shows the correlation between p and the uncorrected multiplicity of tracks 
with 1 7] | < 0.9. The linear increase corresponds to an uncorrected (p t ) of about 0.7 GeV/c per accepted 
charged track. Both (p t ) and multiplicity, and thus also the background p x density, strongly depend on 
the minimum p t threshold (/?™ n ) applied for tracks used as input to the jet finding. To minimize the bias 
on jet fragmentation, a value of p™ m = 0.15GeV/c is preferred. In addition, p™ m = 1 and 2GeV/c are 
investigated to facilitate comparisons to other experiments and to Monte-Carlo generators in a region of 
constant and high tracking efficiency. The mean p over all events and its standard deviation is given for 
different p™ m and two centralities in Table [T] As expected, the mean background p t density decreases 
for larger pf™, for central collisions and p™ m =2 GeV/c it is reduced by an order of magnitude. As one 
can see in the insets of Figure [T] and the standard deviation in Table [TJ the spread of p for the 10% most 
central events is considerable, underlining the importance of the event-by-event background subtraction. 
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Figure 2: 8p t of random cones in the 10% most central Pb-Pb events for the three types of random cone 
probes with p™ m = 0. 15 GeV/c. A Gaussian fit to the left-hand-side and its extrapolation to positive 8p t 
are shown for measured Pb-Pb events and for randomized Pb-Pb events (;it LHS and a LHS in GeV/c). The 
solid line is a fit to the 8 p t distribution for the randomized events with a T distribution shifted to zero 
(Equation [2]) as approximation for the shape in case of independent particle emission. 

4 Sources of Background Fluctuations 

To study the sources of background fluctuations in an unbiased way that is not influenced by a particular 
choice of jet finder, a single rigid cone with radius R = 0.4 is placed in each reconstructed event at 
random <p and tj, with centroid lying within |tj | < 0.5. The background fluctuations are characterized by 
calculating the difference of the scalar sum of all track p t in the cone and the expected background: 

$Pt = Y,Pt,i-A-p, (1) 

i 

where A = 7iR 2 . 

The rigid random cone (RC) area and position are not influenced by the event, so that it provides a 
sampling of the event structure at the typical scale of a jet, but independent of biases induced by the 
choice of a particular jet finding algorithm. The RC measurements will be compared to the embedding 
of jet-like objects, that is directly relevant to the measurement of the inclusive jet spectrum with a specific 
choice of jet finder. 

To characterize the 8p t distribution the standard deviation, a(8p t ), is utilized. In addition, a Gaussian 
distribution with mean jU L and standard deviation cr L is iteratively fit to the distribution within 
[/i LHS - 3a LHS ,/i LHS + 0.5a LHS ], i.e. to the left-hand-side. The a LHS of the fit provides the lower limit 
on the magnitude of the fluctuations and is used to characterize shape differences between the positive 
and negative tails of the distribution, by extr apolating the Gaussian distribution to positive 8p t . 

Various sources contribute to background fluctuations in a heavy-ion event, including: (i) random, un- 
corrected (Poissonian) fluctuations of particle number and momentum; (ii) region-to-region correlated 
variations of the momentum density, induced by detector effects, e.g. a non-uniform efficiency, and by 
the heavy-ion collision itself, e.g. by variation of the eccentricity of the nuclear overlap for collisions 
with finite impact parameter. 

The measured 8p t distribution for random cones in the 10% most central Pb-Pb events is shown in 
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Figure 3: 8p t distribution for random cones, averaged over the full azimuth and separated for three bins 
of random cone azimuthal orientations with respect to the measured event plane. In the bottom panel the 
distributions have been shifted to zero using the mean of the left-hand- side Gaussian fit (jU LHS ). 

Figure [2] The distribution is peaked near zero, illustrating the agreement of the background estimate via 
^-clusters and that due to random sampling of the event. The distribution exhibits an asymmetric shape 
with a tail to the right-hand- side of the distribution, which is also reflected in the difference between 
the standard deviation of the full distribution of a(8p t ) = ll.OGeV/c and the Gaussian width of the 
left-hand-side a LHS (8p t ) = 9.6GeV/c. 

To further differentiate random and correlated sources of fluctuations, three variations of the random 
cone method are employed: (i) sampling of measured Pb-Pb events, (ii) sampling of measured Pb-Pb 
events, but avoiding overlap with the leading jet candidate in the event after background subtraction 
by requiring a distance D = 1.0 in (T],0) between the random cone direction and the jet axis, and (iii) 
sampling of Pb-Pb events in which the (T],0) direction of the tracks has been randomized within the 
acceptance, which destroys all correlations in the event. Figure [2] shows that when avoiding the leading 
jet candidate to suppress upward fluctuations, e.g. due to a hard process, the tail to the right-hand-side is 
already significantly reduced. 

Note that, even for the case of purely statistical fluctuations, the distribution is not expected to be sym- 
metric or to follow a Gaussian shape on the right-hand-side, since the shapes of the underlying single 
particle p t and multiplicity distributions are not Gaussian. In the case of uncorrected particle emission 
a r-distribution provides a more accurate description of the event-wise (p t ) fluctuations [30]. This also 
holds for 8 p t distributions, which are similar to a measurement of (p t ) fluctuations in a limited interval 
of phase space. Taking into account the subtraction of the average background the functional form of the 
probability distribution of 8p t for independent particle emission can be written: 

f(8 Pt ) =A-a h /T(a p ) ■ (a h 8p t + a p ) a "- 1 . e ~^ +a p\ (2) 

This corresponds to a r-distribution with mean shifted from a p /ab to zero and standard deviation a = 
y/ap~/ai,. As seen in Figure |2] this functional form provides a good approximation of the 8p t distribution 
for randomized events, corresponding to uncorrected emission. In this case the distribution is also 
narrower on the left-hand-side. This points to the presence of correlated region-to-region fluctuations in 
addition to purely statistical fluctuations and those expected from hard processes. 

One source of region-to-region variation in the background p t density is the initial anisotropy of the 
nuclear overlap for finite impact parameter collisions, which translates via the collective expansion of 
the medium into an anisotropy in momentum space Q [31] with respect to the symmetry plane of the 
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Figure 4: Dependence of the standard deviation of the 8p t distributions on uncorrected charged particle 
multiplicity, compared to the limit derived from the measured track p t spectrum (Equation [3]> and from 
additional elliptic and triangular flow contributions (Equation |4j). R = 0.4, p™ m = 0.15 GeV/c. 

collision. The event plane direction can be calculated using the azimuthal distribution of all accepted 
tracks within each event, which is dominated by soft particle production. The final state hadron azimuthal 
distribution with respect to the reaction plane of the event, is characterized by a Fourier expansion where 
the leading term is the second moment, called elliptic flow vi. In addition to the geometry driven even 
harmonics (mainly v-i), odd flow components (e.g. V3) driven by intial state fluctuations can modify the 
azimuthal distribution within the event [32]. 

To explore the effect on background fluctuations of azimuthal orientation relative to the reaction plane, 
the 8p t distribution from random cone sampling is studied as a function of the azimuthal orientation 
of the cone axis, (j>, relative to the reconstructed event plane orientation, \j/RP- Three bins are chosen; 
the out-of-plane orientation where the azimuthal angle between the reconstructed event plane and the 
cone axis is > 60°, the in-plane orientation where this angle is < 30°, and the intermediate orientation 
where the angle is between 30 and 60°. The distributions of 8p t for random cones averaged over the 
full azimuth and for the three different orientations are shown in Figure [3] It can be seen that, for out- 
of-plane cones, the most probable background p t density is smaller by almost 5 GeV/c relative to the 
azimuthally averaged estimate of p, with opposite effect in-plane. This shift scales with the average 
flow and the background p t density for a given centrality (oc V2 • p), and is seen to be sizable in central 
events, though discrimination of the event plane orientations is limited by finite event-plane resolution 
ll33l and possible biases due to hard jets. The decreasing width of left-hand-side Gaussian is qualitatively 
consistent with the expectation from reduced particle number fluctuations out-of-plane compared to in- 
plane. For a visual comparison of their shape, the distributions have been shifted such that the centroid 
of the left-hand-side Gaussian fit is zero (see Figure [3]). Notably, the left-hand-side of the distribution 
appears similar for all orientations of the random cones to the event-plane. The random cones distributed 
in-plane show a more pronounced tail to the right-hand-side, compared to out-of-plane. This may point 
to a dependence of the jet spectrum on the orientation relative to the reaction plane, though further 
systematic studies are needed to assess biases in the event plane determination due to jet production 
and possible auto-correlations. For the measurement of the inclusive jet spectrum the correction via an 
event-plane dependent p will reduce the influence of even flow components on the average reconstructed 
jet momentum, but its systematic precision is limited by the finite event-plane resolution. 
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The width of the 8p t distribution due to purely random fluctuations can be estimated from the measured 
single particle p t spectrum via JT8l : 

cr(Spt) = ^N A -a 2 ( Pt )+N A -( Pt ) 2 . (3) 

Here, N A is the expected number of tracks in the cone area A for a given event centrality or multiplicity 
class, {p t } the average p t and o(p t ) the standard deviation of the track p t spectrum. Local variations of 
the average multiplicity, (p t ), or a(pt), lead to additional fluctuations. These region-to-region variations 
can be induced e.g. by (mini-)jets, where the particle p x spectrum is considerably harder than for the 
global event average, and by collective flow. Uncorrelated non-Poissonian (NP) fluctuations can be 
added to Equation [3] knowing their standard deviation, e.g. for additional region-to-region variation of 
the average multiplicity: 

o(S Pt ) = ^N A ■ a 2 { Pt ) + (N A + a 2 P (N A )) • ( Pt ) 2 . (4) 

Figure [4] shows the comparison of the multiplicity dependence of a(8p t ) for the three different types of 
random cones. The distribution of purely statistical fluctuations given by Equation [3] well describes the 
randomized events. Also shown are two parameterizations following Equation[4] Additional multiplicity 
fluctuations due to elliptic flow are approximated from the /^-integrated V2 values measured by ALICE 
for different centralities [7] as a^ P (N A ) 2v 2 N A . This approximate inclusion of V2-effects accounts 
qualitatively for the larger fluctuations in mid central collisions compared to the randomized events 
and the deviation from a v^-increase. The random cone sampling with an anti-bias on the leading 
jet has, by construction, a reduced standard deviation and is close to the parameterization of elliptic 
flow. Taking into account also region-to-region fluctuations from triangular flow, V3, is of particular 
importance in central events where it reaches a similar magnitude as vi IT321 . The contribution of V3 
can be added in quadrature (a^ P (N A ) « 2A^(v| + V3)) since the second and the third harmonic are not 
correlated via a common plane of symmetry [32], for simplicity V3 has been approximated by a constant 
value of V3 = 2.4%. As expected, the inclusion of V3 can account partially for the difference to the 
randomized event in the most central events. In the comparison one has to consider that in practice 
the contribution from hard processes to the right-hand-side tail cannot be cleanly separated from (soft) 
upward multiplicity fluctuations induced by flow. In addition, the approximate description of flow effects 
following Equation|4]does not take into account any flow-correlated changes of (p t ) and o(p t ). 

The track reconstruction efficiency affects the total multiplicity and the shape of the measured p t - 
spectrum at low p t . Using Equation [3] the change of the uncorrelated fluctuations due to finite efficiency 
can be estimated from the efficiency corrected /? t -spectrum in each centrality bin. This procedure sug- 
gests that, for p t mm = 0.15 GeV/c, there is an increase of the standard deviation by 5.4-6.0%, depending 
on centrality. The complete correction requires the knowledge of all correlations within the heavy-ion 
event and is beyond the scope of the present study. 

5 Background Fluctuations in Jet Reconstruction 

The measured jet spectrum in heavy-ion collisions is affected over the entire p t range by background 
fluctuations, especially due to the large and asymmetric tail towards positive 8p t . For the measurement 
of the inclusive jet cross section, background fluctuations can only be corrected on a statistical basis 
via unfolding. Such background fluctuations are evaluated using embedding and reconstruction of a 
probe with identical jet algorithm and parameters as those applied to the data analysis, to account for the 
jet-fmder-specific response to the heavy-ion background. 

In the present study, two probes are embedded into the Pb-Pb events measured by ALICE: (i) single 
high-/? t tracks at various p t , and (ii) pp jet events generated using PYTHIA [34] followed by a detailed 
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simulation of the full detector response. Jet candidates are reconstructed from the event using the anti-& t 
algorithm with R = 0.4 and matched to the embedded probe, by either finding the single track in it, or 
by requiring that the p t of the embedded tracks within the reconstructed jet sum up to at least 50% of the 
original probe jet transverse momentum (p^ mbs ). The difference between the reconstructed, background 
subtracted jet and the embedded probe is then given, similar to Equation [T] by l29l l35ll : 

8p t = p i t et " ec -A^ ec - P - P r be - (5) 

The response may depend on the jet finder, its settings, and the properties of the embedded probe, such as 
p\ m e , area, and fragmentation pattern. In particular the insensitivity to the latter is essential for a robust 
and unbiased reconstruction of jets in heavy-ion collisions, where the fragmentation pattern is potentially 
modified relative to that in pp collisions, and is indeed the observable of interest. 

The dpi distributions measured for each of the methods are shown in Figure [5] Here, the focus is on 
high p^ m e (> 60GeV/c), where the efficiency of matching the embedded probe to the reconstructed jet 
approaches unity. The results are very similar to the random cone method, including the presence of 
an asymmetric tail to the right-hand-side of the distribution. The standard deviations, however, show a 
small increase compared to the random cone method, which is largest for jet embedding (see Table [2]). 
The increase may be due to sensitivity of the jet finder to back-reaction, e.g. the stability of the probe 
area and jet direction after embedding. The single particle embedding can be considered as extreme 
fragmentation leading to rather rigid cones with stable area 7tR 2 , while in the case of true pp-jets the 
probe and reconstructed area may differ, depending on the fragmentation pattern. In addition the finite 
jet area resolution due to the size of the ghost area has to be taken into account [29]. With a ghost 
area of 0.005 a compromise between reasonable jet area resolution and computing time and memory 
consumption was chosen. In the case of track embedding at high p u the jet area resolution fully accounts 
for the difference of 200 MeV/c observed in the standard deviation. 

The broadening of the 8p t distribution for jets with pj™ n = 0. 15GeV/c, as seen in Figure [5] has been 
investigated more closely. The additional left-hand-side structure is caused by probe jets with large area 
(A probe > 0.6) that are split in the heavy-ion event into two separate objects of smaller size. Jets with 
a large area (A > 0.6) are only formed by the anti-& t algorithm in exceptional cases, where there are 
two hard cores at distance close to R Il36ll . It is also seen in Figure [5] that, with increasing p™ m , the 
deviations on the left-hand-side in the case of jet embedding become more pronounced, suggesting that 
the jet-splitting is an effect of hard fragmentation. 

In the determination of 8 p t fluctuations as described above, the probes have been embedded into an event 
population recorded with a minimum-bias trigger. However, the requirement of a hard process biases the 
population towards more central (small impact parameter) collisions, due to nuclear geometry. Correc- 
tion of this effect for centrality bins of 10% width, generates a negligible increase of the fluctuations 
(< 0.1 GeV/c). The full centrality dependence of the fluctuations is given via the standard deviation of 
the distributions and for different p™ m cuts in Table [3] 

The increase of the pj" in cut on the input tracks for jet finding reduces the background fluctuations, due 
to the smaller influence of statistical and soft region-to-region fluctuations. This is observed in Figure [5] 
when the p™ n is varied from 0.15 to 2GeV/c. A p™ m of 2GeV/c reduces the standard deviation by 
more than a factor of two compared to 0.15 GeV/c. Soft region-to-region fluctuations that dominate the 
left-hand- side of the distribution are reduced by a factor of three (see Table [2]). A high pf m significantly 
reduces the impact of fluctuations in the jet spectrum (see Table [4]). However, it may also introduce a 
bias in the jet reconstruction towards hard fragmentation. 

To estimate the influence of the observed fluctuations on the jet measurement, a power law spectrum 
starting at p t = 4 GeV/c has been folded with a Gaussian of width a Gauss = 1 1 GeV/c (5 GeV/c) and 
with the measured 8p t distributions for p™ n = 0.15GeV/c (2GeV/c). The yield increase relative to the 



Background Fluctuations for Jet Reconstruction in Pb-Pb collisions 



9 



| 1 

| 10- 1 
|mo- 2 

1 10- 4 
°-w 5 



1 

Iff 1 
10" 2 
10" 3 
10" 4 
10" 5 



1 

10" 1 
10" 2 
10 3 
10" 4 
10" 5 



-50 50 100 

6p ( ch (GeV/c) 

Figure 5: 8p t distribution of charged particles for jet reconstruction with the three methods in the 10% 
most central Pb-Pb events for pf n = 0.15 GeV/c, 1 GeV/c, and 2 GeV/c. 

unsmeared spectrum in one high ;? t -bin for the most central collisions is given in Table|4] For the different 
probes, they agree within the uncertainties given by the statistical fluctuations in the tails of the 8p t - 
distributions; about a factor of ten increase for p™ m = 0.15 GeV/c and a 30% effect for p™ n = 2GeV/c. 
Minor differences in the standard deviation as well as the left-hand-side differences have no sizable 
effect on the spectral shape after folding. The difference between smearing with the full 8 p t and with 
a Gaussian distribution illustrates the strong influence of the right-hand-side tail, which must be taken 
into account in the analysis of background fluctuation effects on jet reconstruction. The extracted values 
naturally depend on the choice of the input spectrum, so, in addition to the power law, a jet spectrum 
for pp collisions at ^/^nn = 2.76 TeV extracted from PYTHIA simulations has been used. These studies 
indicate that the increase of yield due to background fluctuations falls below 50% for reconstructed 
charged jets in the region of p t ^ 100 ± 15GeV/c (60 ± lOGeV/c) for pf n = 0.15GeV/c (2GeV/c). 
Repeating the exercise with a Gaussian smearing of a a Gauss = 1 1 GeV/c and with the 8 p t distribution of 
random cones avoiding the leading jet for p™ m = 0.15 GeV/c leads, as expected, to a reduced influence 
of the tail. Here, the relative yield increase falls below 50% in the range of p t « 75 ± lOGeV/c. The 
employed input spectra do not consider the geometrical limitation on the number of jets that can be 
reconstructed within the acceptance for a single event OTTl . This effect also limits the extraction of jet 
spectra at lower p t via unfolding. 

6 Summary 

The first detailed study of event background fluctuations for jet reconstruction using charged particles 
in Pb-Pb collisions at the LHC has been presented. The standard deviation of the fluctuations in the 
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a (GeV/c) 


a LHS (GeV/c) 


H LHS (GeV/c) 


A min = 0.15GeV/c 


random cones 


10.98±0.01 


9.65 ±0.02 


-0.04 ±0.03 


track emb. 


11.19±0.01 


9.80±0.02 


0.00 ±0.03 


jet emb. 


11. 34 ±0.02 


9.93 ±0.06 


0.06 ±0.09 


pf B = 1 GeV/c 


random cones 


8.50±0.01 


7.08±0.01 


-0.22 ±0.02 


track emb. 


8.61 ±0.01 


7.11 ±0.01 


-0.25 ±0.02 


jet emb. 


8.78 ±0.02 


7.25 ±0.04 


-0.08 ±0.08 


pf B = 2 GeV/c 


random cones 


4.82±0.01 


3.41 ±0.01 


-0.01 ±0.01 


track emb. 


4.88±0.01 


3.05±0.01 


-0.92 ±0.01 


jet emb. 


5.03 ±0.01 


3.52±0.01 


0.01 ±0.02 



Table 2: Background fluctuations in central events. Comparison of the Gaussian fit to the left-hand- 
side of the S /^-distributions and the standard deviation in central Pb-Pb collisions for the three different 
methods and for the three /?™ n -cuts. The quoted uncertainties are purely statistical. 







pf n (GeV/c) 






0.15 


1.0 


2.0 


Centrality Class 


o 


(5a) (GeV/c) 




0-10% 


11.19±0.01 


8.61 ±0.01 


4.88±0.01 


10-20% 


10.19±0.01 


7.67±0.01 


4.29 ±0.01 


20-30% 


8.46 ±0.01 


6.35±0.01 


3.58±0.01 


30-40% 


6.51 ±0.01 


4.93±0.01 


2.68±0.01 


40-50% 


4.71 ±0.01 


3.63±0.01 


1.95±0.01 


50-60% 


3.28±0.01 


2.61 ±0.01 


1.41 ±0.01 


60-70% 


2.22 ±0.01 


1.70±0.01 


0.95 ±0.01 


70-80% 


1.48 ±0.01 


1.01 ±0.01 


0.62±0.01 



Table 3: Centrality dependence of fluctuations. Standard deviation of 8p t distributions and statistical 
uncertainty for different centrality bins and p™ m cuts using the track embedding probe. 



/(a) folded with 


relative yield for p t = 60 


-68GeV/c 


8p t 


RC 


tracks 


jets 


p™ m = 0A5GeV/c 


9.8±1.7 


11.4± 1.1 


10.9 ±3.4 


pf B =2GeV/c 


1.30 ±0.02 


1.31 ±0.02 


1.65 ±0.25 


Gauss 








a= HGeV/c 




1.82 ±0.04 




(7 = 5 GeV/c 




1.05 ±0.01 





Table 4: Yield modification for power law spectrum. Relative yield in the bin p t = 60 — 68 GeV/c for a 
power law spectrum (f{pt) = 0.7/(0.7 + and p t > 4GeV/c), folded with the different Sp t distribu- 
tions for 0-10% centrality and with a Gaussian, where the width is similar to the standard deviation of 
the Sp t distributions. 
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10% most central events is a = (10.98 ±0.01) GeV/c within a rigid cone of R = 0.4 and for a low 
p t cut-off of 0.15GeV/c. It has been shown that the non-statistical sources of fluctuations are driven 
in part by the anisotropy of the particles emitted from the collision (elliptic and triangular flow). The 
variation of multiplicity in different orientations with respect to the event plane, induces shifts in the 
background-subtracted jet p t even for central Pb-Pb-collisions. 

The anti-& t jet finder response for charged particle jet reconstruction has a modest dependence on the 
method used to characterize the fluctuations. For embedded, simulated pp-jets the standard deviation 
increases to (1 1.34 ± 0.02) GeV/c. In addition, certain rare fragmentation patterns in pp are likely to be 
split in the heavy-ion environment leading to minor effects in the background response. The observed 
differences between the two types of embedded probes (namely single tracks and pp jets) do not indicate 
a strong sensitivity of the reconstructed anti-& t jet spectrum on fragmentation. The case of a strong 
broadening of the jet due to medium effects has not been considered here. 

The use of reconstructed charged particles down to p™ m = 0.15GeV/c allows a comparison of the 
impact of background fluctuations with a minimal bias on hard fragmentation in jet finding to the 
case with increased bias (p™ m > 1 GeV/c). The observed reduction of the standard deviation to a = 
(4.82 ± 0.01) GeV/c for the unbiased sampling and p™ m = 2GeV/c is driven by the smaller number 
fluctuations and the reduced influence of soft region-to-region fluctuations. 

The asymmetric shape of the 8p t distribution with a tail towards positive fluctuations has a large im- 
pact on the jet measurement, compared to purely Gaussian fluctuations, though the role of signal jets 
contributing to the tail has to be considered. Using different assumptions on the shape of the true jet 
spectrum it is found that for p™ m = 0. 15 GeV/c fluctuations can have a large influence on the charged jet 
yield for transverse momenta up to 100 ± 15 GeV/c. 
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